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Abstract: The recent experimental realisation of a one-dimensional Bose gas 
of ultra cold alkali atoms has renewed attention on the theoretical properties of 
the impenetrable Bose gas. Of primary concern is the ground state occupation 
of effective single particle states in the finite system, and thus the tendency for 
Bose-Einstein condensation. This requires the computation of the density matrix. 
For the impenetrable Bose gas on a circle we evaluate the density matrix in terms 
of a particular Painleve VI transcendent in cr-form, and furthermore show that 
the density matrix satisfies a recurrence relation in the number of particles. For 
the impenetrable Bose gas in a harmonic trap, and with Dirichlet or Neumann 
boundary conditions, we give a determinant form for the density matrix, a form 
as an average over the eigenvalues of an ensemble of random matrices, and in 
special cases an evaluation in terms of a transcendent related to Painleve V 
and VI. We discuss how our results can be used to compute the ground state 
occupations. 



1. Introduction 

Recent advances in the experimental physics of Bose-Einstein condensates |14l 
I15ll^ have led to the experimental realisation of a one-dimensional Bose gas of 
ultra-cold alkali atoms. One expects [211 that the microscopic forces are such 
that there is an effective one-body confining harmonic potential acting on each 
atom individually, and an effective infinitely short range contact potential acting 
between neighbouring atoms. Moreover, in a certain physical regime depending 
on the ratio of the transverse confinement width to the s-wave scattering length, 
it is argued in '39' that the contact potential can be well approximated by the 
delta function form U{\x — y\) = gS(\x — y\), and furthermore (7 — > oo in the low 
energy scattering limit. The limit 5 ^ 00 of the delta function interaction Bose 
gas is the impenetrable Bose gas, introduced in [T ^ .p ^ . 
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Not surprisingly, there has thus been renewed interest in the theoretical prop- 
erties ol the ground state of the finite system impenetrable Bose gas Kffllll3| . With 
the 3d Bose gas exhibiting Bosc-Einstein condensation, a central question is the 
tendency of the finite system confined to Id to form a Bose-Einstein condensate. 
To attack this question is a two step process. First, with the particles confined 
to the region i? € M and the ground state wave function tpo real, it is necessary 
to compute the one-particle density matrix 



PN{x-y)=N dx2.-.j dxNi^oix,X2, ■ ■ ■ ,XN)i^o{y,X2, ■ ■ ■ ,xn), (1-1) 
Second, one must solve the eigenvalue problem 

PN{x;y)(t>k{y)dy = Xk(t>kix), /c G Z>o. (1.2) 



f2 

Because this integral operator is idempotent, the Xj are non- negative, while 
the trace condition J^^ pn{x;x) dx — N implies J^k'^k — N. Consequently the 
Afc have the interpretation as occupation numbers of effective single particle 
states (j)kix). The simplest case is when ^2 = [0, L] with periodic boundary 
conditions. The periodicity implies that pn{x] y) — pn{x — y; 0). Thus we have 
0fe(x) = -^e^'^^'^^^Z-f^ and so 

Xk= [ pAr(a;;0)e2-*'="/^dx. (1.3) 







However for other geometries and confinements there is no analogue of 1)1. 3|l and 
one must solve (|1.2|l numerically. 

A number of results are available on p]y{x;0) for periodic boundary con- 
ditions. In particular Lenard has given pj^^i{x;0) as an N x N Tocplitz 
determinant (see 1)2.511 - 1)2. 15|l below), and subsequently obtained the N ^ oo 
asymptotic expansion 



p^{x;0)^PoA[——^ — ) , A = ^!|/^ (1.4) 



1/2 ^ ^ ^ 

Nsm{npox/N) J ' ^ 

where po denotes the bulk density and G{x) denotes the Bairn's G-function, valid 
for x/N fixed. Although the analysis of leading to 1)1. 4|l was not rigorous, the 
setting of the problem as belonging to the asymptotics of Toeplitz determinants 
with symbols having zeros on [0, 2tt) was identified, and this work inspired a 
subsequent rigorous proof |46j . (We remark that the asymptotic form of Toeplitz 
determinants of this type was first conjectured by Fisher and Hartwig [7], ) 
The result (|1.4|) substituted into 1)1. 3(1 gives Aq ~ cy/N for a specific c computable 
from 1)1. 4|) . Thus for large N the fraction of particles in the zero momentum 
state is proportional to \/N. The result ()1.4|) can also be used to compute the 
large TV behaviour of Afc for any fixed fc > For the impenetrable Bose 
gas confined by a harmonic one-body potential, or indeed in other geometries 
such as Dirichlet or Neumann boundary conditions, no results of this type are 
known. All one has is the recent numerical study of Girardeau et al. in the 
case of the harmonic well, who by a Monte Carlo study of system sizes up to 
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= 10 obtained the estimate Aq cx iV°-^^ for large N. If correct, this resuh 
imphes the maximum effective single particle state occupation is dependent on 
the geometry /confining potential. 

To further study this issue, we take up the first step in the procedure above 
to compute the Xj, and thus provide formulas suitable for the numerical com- 
putation of pN{x\y). Four cases are considered — when the domain is a circle 
(or equivalently periodic boundary conditions) ; a line with the particles confined 
by a harmonic one-body potential; and an interval with Dirichlet or Neumann 
boundary conditions. The Toeplitz determinant formulation in the case of peri- 
odic boundary conditions is extended to Hankel determinant forms for pN{x]y) 
in the other cases (Section 2.2), and a formulation for efficient Monte Carlo 
evaluations by way of expressing the pN{x;y) as averages over the eigenvalue 
probability density function (p.d.f.) of certain matrix ensembles is given (Sec- 
tion 2.3). We then give a systematic Fredholm type expansion of Pn{x\ y) about 
the density p]y{x;x) (Section 2.4). 

Beginning in Section 3 we address the issue of closed form evaluations of 
Pn{x] y). In the infinite system there are some celebrated instances of such eval- 
uations. In particular Jimbo et al. |22j related the problem of evaluating Pooix; 0) 
to integrable systems theory, and consequently were able to derive the formula 

Poo(a;;0) = poexp (^y cry(t)yj, (1.5) 

where ay satisfies the non-linear equation 

{xa'^f + A{xa'y - ay - 1) [xa'y - ay + {a'yf'^ = (1.6) 
subject to the x ^ boundary condition 

av{x)^^^-^ + ^ + 0{x*). (1.7) 

The differential equation l|1.6|) is an example of the so-called Jimbo-Miwa-Okamoto 
(T-form of the Painleve V equation, the latter being essentially the differential 
equation obeyed by the Hamiltonian in the Hamiltonian formulation of Py |38| . 

4 

{th'{rf - [hy - th'y + 2{h'y ff + ^llih'y + Vk) ^ (1.8) 

k=l 

with wi + t;2 + U3 + f4 = 0. Setting 

<Jv(x) + y2^hv{t), = (1.9) 

shows that (|1.6|) reduces to H1.8|l with {vi,V2, ^3,1)4) = (V2, — V2, ^2^—^/2)- Subse- 
quently the characterisation of Paoix] 0) in terms of the solution of a differential 
equation was extended by Its et al. 1211 (see also |2H|) to the characterisation of 
p^(x\0) — the density matrix of the impenetrable Bose gas at non-zero tem- 
perature T, as the solution of coupled partial differential equations. 

In the same study that (|1.6|l was obtained, Jimbo et al. evaluated the scaled 
probability of an eigenvalue free interval for large GUE random matrices (random 
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Hermitian matrices) in terms of another particular case of the cr-form of Py ■ In 
recent years there has been considerable progress in the evaluation of probabil- 
ities and averages in matrix ensembles in terms of Painleve transcendents (see 
e.g. Because of the close relationship between the density matrix for im- 

penetrable bosons and gap probabilities in matrix ensembles, the random matrix 
results can be used to extend the density matrix Painleve transcendent evalu- 
ation (|1.5(l to the exact Painleve transcendent of pn{l{x);x) in the four cases, 
where l{x) denotes the image of x reflected about the centre of the system. 

We adopt two distinct strategies to obtain the exact evaluations. In Section 
3 we present the first approach where we work directly with the definition of 
Pn{x; y) on a circle as a multidimensional integral. It turns out that this integral 
is one of a general class which have recently ^Jj been identified as r-functions 
for certain Pyi systems. We show that our Pyi transcendent evaluation for the 
finite system scales to the infinite system result (jl.6(l . As well as being a special 
case of the class of integrals related to Pyi systems in 11 , the multidimensional 
integral formula for pn{x] y) on a circle is also a special case of a class of integrals 
over the unitary group shown to satisfy integrable recurrence relations in ,1 . We 
will show that these recurrences can alternatively be derived from orthogonal 
polynomial theory ^3,3,- 

Underpinning the second of our strategies is the formulation of Lenard |29j 
which allows pn{x] y) to be expressed in terms of the Fredholm minor of 1 — ^Kj, 
where Kj is the integral operator on J = [x, y] with kernel K of Christoffel- 
Darboux type. It is this formulation which also underlies the calculation of |22| . 
The Fredholm minor in turn can be expressed in terms of the product of the 
corresponding Fredholm determinant, and the resolvent kernel R{s,t) evaluated 
at the endpoints x,y of J. These latter two quantities have been extensively 
studied in the context of gap probabilities in random matrix ensembles ^42 . 43"|l49l 
I48| . allowing us to essentially read off from the existing literature an expression 
for pff(i{x);x) in terms of Painleve transcendents in each case. This is done in 
Section 4. 

The significance of our results, from the viewpoint of the theory of the ground 
state occupation of single particle states for the impenetrable Bose gas, and from 
the viewpoint of the Painleve theory, is discussed in Section 5. 

2. Formulations of pN{x;y) 

2.1. The wave functions. We will first revise the construction of the ground state 
wave function for impenetrable bosons on the circle, on the line with a confining 
harmonic potential, and on an interval with Dirichlet or Neumann boundary 
conditions. The wave function and density matrix will be given a superscript 
"C", "H", "D" and "N" respectively to distinguish the four cases. 

In general the wave function ■0(2;i, . . . , s^at) for impenetrable bosons must 
vanish at coincident coordinates. 



and satisfy the free particle Schrodinger equation otherwise. But for point parti- 
cles without spin the condition (|2.1(l is equivalent to the Pauli exclusion principle. 
This means that for any fixed ordering of the particles, 



V'(a;i, ...,Xi,... 



,Xj, . . . ,xn) = foTXi = Xj, {i ^ j), 



(2.1) 



Xi < X2 < ■ . ■ < Xn 



(2.2) 
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/ dxi ■ ■ 


• / dxN [ 


Jo 


Jo 



say, there is no distinction between impenetrable bosons and free fermions jl2| . 
Consequently the ground state wave function ipo can, for the ordering H2.2|l . be 
constructed out of a Slater determinant of distinct single particle states. For 
other orderings i/jq is constructed from the functional form for the sector (|2.2|l 
by the requirement that it be a symmetric function of the coordinates. 

Consider the case that the particles are confined to a circle of circumference 
length L. This means we require 

^Pixl, ...,Xi + L,.. .,xn) = -^{xi, ...,x.i,... ,xn) (2.3) 

for each i = 1, ... ,7V. Constructing a Slater determinant obeying H2.3|l out of 
distinct single particle states with zero total momentum and minimum total 
energy gives 

, fdetfe^'^^'^^"^/^] ,=i,....« iVodd 

„/,C/„ ^ ^ _ ( i\T\\ — 12 T —N 2 ) fc=-(jv-i)/2,...,(jv-i)/2 

V^o(a.i,...,x^)-(iV!) 'L M 2..(fe+i/2).,A] Neven 

\ k = -N/2,-.., N/2-1 

= (7V!)-i/2l-^/2 Y[ 2ismTr{xk-Xj)/L (2.4) 

l<j<k<N 

where the factor of (iV!)^^/^ is included so that 

b^ixi, . . .,xn)\^ = 1. 

Excluding the (unitary) factors of i, and recalling H2.2(l . we note that this state is 
non-negative — a property which distinguishes the ground state in Bose systems. 
By the requirement that the wave function for a Bose system be symmetrical 
with respect to interchanges Xj <-*■ Xj' {j ^ j') we see immediately from (|2.4|l 
that for general ordering of particles 

V'^(xi,...,a;Ar) =L-^/2(iV!)-i/2 [] 2| sin^(xfe - x,)/L|. (2.5) 

l<j<k<N 

In the case of impenetrable bosons on a line with a confining harmonic po- 
tential, we take as the Schrodinger operator (in reduced units) 

-E^+E4 (2-6) 

The corresponding normalised single particle eigenstates {4'k{.x)\k=o,i,... have 
the explicit form 

0fc(x) = ^e-^'/2^,(x), (cH)2 = 7ri/22-^-fc! (2.7) 

where Hk{x) denotes the Hermite polynomial of degree k. Forming a Slater 
determinant from the minimal energy states (/c = 0, 1, . . . , iV — 1), making use 
of the Vandermonde determinant formula 

det[pj^i{xk)]j.k=i,...,N = det[xlr^].i,k=i,...,N ^ Y]_ {xk - Xj) (2.8) 

l<j<fe<JV 
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for any {pj{x)} with Pj{x) a monic polynomial of degree j, and arguing as in 
going from (j2.4(l to H2.5|l shows 



1 ^ 

N j=i l<j<k<N 



N-1 



[C^f = N\J{ {off 



1=0 



(2.9) 

Finally we consider the case of impenetrable bosons on the interval [0, L] with 
Dirichlet or Neumann boundary conditions, requiring that the wave function or 
its derivative vanishes at x = 0, L respectively. The single particle eigenstates 
{0fe(a;)} in these situations are, in increasing order of energy. 



^sin^, (fc = l,2,...), 



1 

I' 

2 Trkx 



Recalling the C and D type Vandermonde formulas 001 



det[z^^ - %,k= 



life -2/) n ^^^"^^H^^i^J 

j=l l<j<k<n ^ ^ 

--l,....n = 2 TT (Zfe - Zj) f 1 - — 



l<j<k<n ^ ■' 

we see that the corresponding ground state wave functions are 



'tp^{xu...,XN) = — L=| det[(/)^(a;j)]j,fc=i,...,Ar 



1 / 1 



Nl \V2L 



N N 



Jl2sin(7rx//-L) 2| cos7rxfc/L — cos7ra;j/i|, (2.10) 

l<j<k<N 



1=1 



and 



2po{xi,...,XN) 



1 



1 1/1 



det[(l>k_i{xj)]j^k=i,...,N 

N-l 



H 2\ COSTTXk/ L — COSTTXj / L\. (2.11) 

l<j<*;<Af 



2.2. The density matrix as a determinant. The density matrix pn+i is defined 
as an A''-dimensional integral by In the cases of the impenetrable Bose 

gas wave functions of the previous section, this integral can be reduced to a 
computationally simpler iV-dimensional determinant. For the circular case, this 
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form has already been given byLenard Thus using the general Heine identity 



TVIdet 



1/2 



dx w{z)z^ ^ 



-1/2 



.1/2 .i/z - 

/ dxi---l dxNWw{zi) W 

J -1/2 J -1/2 .^tt. 



J j.,k=l,...,N 
»l/2 N 



^1/2 1=1 l<j<k<N 

we see from H2.5|l and that 

PN+ii^:^) = 2det[af_k{x)]j^k=i,...,N 

,1/2 



\z, - z, := e^-"^/^ 



J-l/2 



Furthermore, the elements af have the explicit evaluation 



(2.12) 



(2.13) 
(2.14) 



C _ ± ±r^x/L 



TTX TT 2x, TTX 



, 2x . 2ttx 
,(l__)+sm — 



,C _ * \ ^) ±irmrx/L 



TT m{m? — 1) 



Tree . m'KX . TTX mnx 

cos — sm — m sm — cos — - — 



, |m| > 1 
(2.15) 



In particular, it follows that 
1 



Pi(a^;0) = 



X , ,Trx , 



,TTX , 



P2(a;;0) = _[7r(V2- -)cos( — ) +sin( — )] 
8 



(2.16) 
(2.17) 



p^ix; 0) = ^{2 - y.nHy^ + MV2 - f ) sin(^) cos(^) 

+ h% + 2n^y, ~ cos2(^) + V2Cos4(^)}. (2.18) 

An essential ingredient underlying the applicability of (|2.12|) is the factorisa- 
tion 



(a;,xi, . . . ,xjv) 



1 1 



AT 



observed from (|2.5|l , and their subsequent use of the determinant form in 
to replace tAq' on the right hand side. Now we observe from (|2.9I) . H2.10|l . H2.11|l 
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that -0^^, ip^ , ipQ in the case of + 1 particles can similarly be factorised. Using 
the general identity 



mdet 



dt g{t)h,^i{t)hk-i{t) 

f°° f°° -fr 2 

= / dxi--- dxN [[gixi) {det[hj_i{xk)]j,k=i,...,NY 



j,k=l,...,N 
N 



1=1 



(c.f. (|2.12|l ') we thus obtain analogous to (|2.13I2.14|I the determinant formulae 



Pn 



Pn+A^'^ y)^2 T T °*^^K\fc(^! y)\]M=i,...,N 

PN+i{x;y) = ■^det[a2f^{x;y)]j,k=i....,N, {N > 1) 



where 



2— j — fc+2 pco 



(2.19) 

(2.20) 
(2.21) 

(2.22) 



a^i^{x;y)^8 / | cos cos7rt||cos cos 7rt| sin Trjt sin Trfci (2.23) 

Jo L L 

I TTX '^y 

a^i.i^;y)—8 / | cos — cos7ri||cos— cos7rt| cos7r(j — l)i cos7r(/c— l)t. 

Jo L L 

(2.24) 



To simplify further, we note 
\x - t\\y - t\ 



{x-t){y~t), t^[x,y] 
-{x - t){y - i), te [x,y\, 



(2.25) 



and similarly with | cos ^ — cos7ri||cos^ — cos7rt|. Use of such an identity 
allows aj^^(x; J/) to be evaluated in terms of incomplete gamma functions, and 
a}^j.{x\y), a^i,{x;y) in a form similar to ()2.15f) . 



2.3. pM+i{x]y) and integrals over the classical groups. In general, for a many 
body wave function ■f/^o, iV'oP has the interpretation as a multivariable p.d.f. As 
first observed by Sutherland in the cases of V'o^ and -^c?, a feature of It/^oP for each 
of the wavefunctions H2.5|l . (|2.9() . H2.10|l and p. 11(1 is that it coincides precisely 
with the multivariate p.d.f. for particular classes of random matrices. Thus 

|Vo'^P=Ev(C/(7V))|,^2Wi 
|0?P =Ev(GUEAr) 

|Vo°l'=Ev(5p(iV))|,^,,/i 

100^12= Ev(0+(2iV))|,=.,/i 
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where Ev{X) denotes the eigenvalue p.d.f. of the ensemble of matrices X, and 
U{N) denotes the unitary group with uniform (Haar) measure, GUEjv the Gaus- 
sian unitary ensemble of NxN complex Herniitian matrices, Sp{N) the group of 
symplectic unitary 2N x 27V matrices with Haar measure, and 0+ (2 A^) denotes 
the group of real orthogonal 2N x 2N matrices with determinant +1 and Haar 
measure. 

Moreover, it follows from the definition Hl.l(l of the density matrix, and the 
explicit forms of the wave functions, that pN+i{x; y) in each of the cases can be 
written as an average over Ev(X) for appropriate X. Explicitly 



P%+ii^;0) - 7 (n - |)l|2sin(|)|^ (2.26) 



Ev(C/(Ar)) 
N \ 



Ev(GUEjv) 



PN+i(a;;2/) 



2 Txx Try 
— sin — sin — 
ij Li ij 



X / J|2|cos^ -cos6';|2|cos^ - cose*,! \ (2.28) 

\'=1 / Ev(Sp(iV)) 

P^^r+i(a;;2/) = ^/n2|cos^-cos0i|2|cos^-cos0/|\ . (2.29) 

\'=1 / Ev(0+(2Ar)) 

Because it is straightforward to generate typical members from each of these ma- 
trix ensembles (see e.g. |^), and so compute eigenvalues from the corresponding 
p.d.f. Ev(X), these expressions are well suited to evaluation via the Monte Carlo 
method. 

For future reference we note that the density matrices for the corresponding 
free Fermi systems are given by the same averages, except that the absolute 
value signs are to be removed. In particular 

p^^^f (x;0) = \ (n2sin(| - f )2sin(|)\ . (2.30) 

Furthermore it is elementary to compute density matrices for free Fermi systems 
as sums over single particle states (a consequence of all energy states below the 
Fermi surface having occupation unity), and this implies the explicit evaluation 

C,FF, 1 sin(7r(jV + l)x/L) 

"^+i("'°)=L sin(W^) • ^ ^ ^ 



2.^. Systematic small-\x — y\ expansion of pN{x]y). According to the definition 
the density matrix at coincident points x = y is equal to the particle 
density. But the particle density for the impenetrable Bose gas is the same as for 
the corresponding free Fermi system and thus simple to compute. In the infinite 
system, the translational invariance of the state gives that the particle density 
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is a constant. For this case Lenard 30 has shown how to make a systematic 
expansion of the density matrix pao{x;y) about the case of coincident points 
Pooix'jx). Here we will present this expansion for finite Bose gas systems with 
ground state wave functions of the form 



1 ^ 

~n^(^') n \xk-xj\. 



c 



(2.32) 



1=1 



l<j<k<N 



This form includes the case of the harmonic well l|2.9|l . and after the change of 
variables cos irxj /L^ Xj in H2.10|l and H2.11|l also includes the case of Dirichlet 
and Neumann boundary conditions. 

Following Lenard ISO], we note that substituting H2.32|l in and using 

shows 



, ^ Ng{x)g{y) ( 
PN{x;y) = 



dx2 g'^{x2) ■ 



dxN — ^ 



N 



Y[{x - xi){y - xi) Y[ {xk-Xj f 

2<j<k<N 



1=2 



?=2 



dxN g'^{xN) 
(2.33) 



One now introduces the Fermi type distribution function 
Ng{x)g{y) 



p^{x;y;x2, ■■■,Xn) 



C 



■Y\_9^(^l) / dXn+l g'^{Xn+l) 
1=2 •' ^ 



N 



Y\{x - xi){y - xi) Y\ {xk~XjY 

2<j<k<N 



dxN g^{xN) 
(2.34) 



1=2 



(when n = 1 this corresponds to the free fermion one-body density matrix). 
Expanding H2.33|) in a power series in ^ and using the definition 12.34|l shows 



PN{x;y) = Y^ — dx2--- 

n=o "'^ "'^^ 



dXn+lPN [X] y]X2,.--, Xn+l) 



5=2 



(2.35) 



(the summation can be extended to infinity since p^^ {x;y;X2, ■ ■ ■ ,Xn) = for 
n>N). 

Next, let {pj(a;)}j=o,i,... be monic polynomials of degree j, orthogonal with 
respect to the weight function g'^{x). Then writing the integrand in (|2.34(l as a 
product of Slater determinants using (|2.8(l and making use of the orthogonality 
of the pj (x) , a standard calculation shows 



pfl{x;y;x2, ... ,Xn) = det 



=: K 



K{x,y) [K{xj,y)]j=2,...,n 

[K{x, Xk)]k=2,...,n [K{Xj,Xk)]j,k=2,...,r 
X ' ' ' X j-j 

,y X2 ■■■ x„ 



(2.36) 
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where, with 7V^ := g'^{x){pj{x))^ dx 



K{x,y) := g{x)g{y) ^ 



9{x)g{y) pn{x)pn-i{v) - PN~i{x)pN iy) 



(2.37) 



Mv-i x-y 

The equahty in (|2.37|) follows from the Christoffel-Darboux summation formula, 
and leads to the name Christoffel-Darboux kernel (the latter term is due to a 
relationship with integral equations; see Section 4.1) for H2.37|l . Hence 



Pn 



{x; y) = ^ — / dx2 - ■ ■ dxn+i K 

n=0 "'^ "'^ 



X X2 ■ ■■ Xn+1 
y X2 ■■■ Xn+1 



C=2 



i=2 



(2.38) 



where 



n=0 •'^ 

As for — y| small each term in (|2.38ll is proportional to successively higher 
powers of |a; — ?/|, this is the sought systematic small |a; — y| expansion of Pn{x] y)- 
We will see in Section 4.1 that the expansion (|2.38|) forms the basis for Painleve 
transcendent evaluations of pn{i^{x), x) in the harmonic well, Dirichlet and Neu- 
mann boundary condition cases. 



3. Jimbo-Miwa-Okamoto r-functions and orthogonal polynomials 

In this section we will provide the finite N analogue of the Jimbo, Miwa, Mori and 
Sato |22) Painleve transcendent evaluation H1.5|l of Poc{x] 0), by similarly evalu- 
ating p%j^i{x\ 0), and also presenting a recurrence relation in N for p^j^-^{x] 0). 
Our Painleve transcendent evaluation of p^j^^{x]0) is in terms of the solution 
of the Painleve VI equation in a-form. Let us then discuss some of the theory 
relating to this equation. 



3.1. Hamiltonian formulation of Pyi and r-function sequences. There are six 
Painleve equations, labelled Pj - Pyi . They result (see e.g. 18 ) from the project 
undertaken by Painleve, Gambler and others to classify solutions to second order 
differential equations of the form y" = R{y',y^t), where R is rational in y', 
algebraic in y and analytic in t which are free from movable branch points. It 
was shown that the only such equations, excluding those which could be reduced 
to first order equations or to linear second order equations, are Pj - Pyi . Our 
interest is in the Pyi equation, which has the form 

\q q-1 q-tj \t t-1 q-tj 

q{q-l){q~t)f ^ {t - I) , , t{t - l) \ 

'a + /?^+77 Tv2+0l IT? h (3.1) 



^2(^-1)2 ^q2 '(g-l)2 {q-ty 
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and its solution, the Pyi transcendent q{t). We will see that p'^^^{x;0) can be 
identified with a r-function sequence in the Pyi system. 

The Pyi system refers to the Hamiltonian system {q,p; H, i\ 

q' -^^ p' - (3 2) 

dp ' dq 

where, with ao + Q^i + 2q!2 + 0:3 + 0:4 = 1, 

t{t-l)H = q{q-l){q-t)p' 

- [ai{q - l){q - t) + asqiq - t) + (ao - l)q{q - l)]p 

+ a2{ai + a2){q - t). (3.3) 

It has been known since the work of Malmquist in the early 1920's 34 that the 
Pyi equation l|3.1|l results from the Hamiltonian system 1)3. 2|l . 1)3.3(1 by eliminat- 
ing p and choosing the parameters so that 

a = Vaa?, /? = -^^al 7 = V2al, S = 1/2(1 - al). 

One sees that the Hamiltonian can be written as an explicit rational function 
of the Pyi transcendent and its derivative. This follows from the fact that with 
H given by H3.3I) . the first of the Hamilton equations is linear in p, so p can be 
written as a rational function of q, q' and t. 

The r-function is defined in terms of the Hamiltonian by 

if=|logr(t). (3.4) 

The utility of being able to identify p' ^_f_^{x ; 0) as a r-function for the Pyi system 
is that H, and thus by integration of 1)3. 4|l r(i), can be characterised in terms of 
a differential equation. 

Proposition 1. \21[{3T\I Rewrite the parameters ao, . . . , 04 of iS.'-^) in favour of 
the parameters 

bi = y2(a3+a4), 62 = V2{a4-a3), 63 = y2(ao+cki-l), 64 = ^^(ao-ai-l), (3.5) 
and introduce the auxiliary Hamiltonian h by 

h:^t{t- l)H + e'2[h]t - ^e2[b] 

= i(t -l)if+ (6163 + 6164 + M4)i-^ ^i^k, (3.6) 

l<J</£<4 

where e^ [b] denotes the jth degree elementary symmetric function in 61,63 and 
64 while ej[h\ denotes the jth degree elementary symmetric function m 61, . . . , 64. 
The auxiliary Hamiltonian satisfies the Jimbo-Miwa-Okamoto a-form of Pyi 

h'(t{l - t)h!'^ ^ + {]i!\2h - {2t - l)h'] + 61626364) ^ = Y[{h' + bl). (3.7) 

k=l 
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A self contained derivation of this result can be found in JI] . 

One of the main practical consequences of the Hamiltonian formulation is that 
it allows for a systematic construction of special solutions via Backlund trans- 
formations - birational mappings which leave the Hamilton equations formally 
unchanged |37j . The elementary Backlund transformations form an extended 
affine Weyl group of type D^P . By composing certain of these elementary opera- 
tors, shift operators can be constructed which have the effect of incrementing the 
a parameters by ±1 or 0. For example, one such operator of this type, denoted 
T3 in Jl], has the action 

T3Q; = (ao + l,ai -I- l,a2 - l,a3,Q;4) 

or equivalently, after recalling (|3.5I) 

Tab = (61,62,63-1-1,64). (3.8) 

Although acting on p and 5 is a non-trivial rational mapping, when acting on 
H , has the formal action of acting only on the a's, 

TiH = H 

This motivates introducing a sequence of Hamiltonians 

T^H ^ H 

and a corresponding sequence of r-functions specified by 

Tg"// = ^lograM, r3[n] =r3[n](i) = r(i; 61, 62, 63 + ^,64). (3.9) 

A crucial result due to Okamoto [23, which can be derived from the specific 
form of the action of T3 and Tg"^ on H, p and q |2Sj, is that r^ln] satisfies a 
particular differential recurrence relation. 

Proposition 2. The r-function sequence 1^3. 9\) satisfies the Toda lattice equation 
^^logr3[n] = M!iZ^p[!i±Il, S = tit-l)i (3.10) 

where 



To Inl ' dt 



/ s (n+&i+b3)(n-(-f)3-(-&4)/2 

f3[n] := (t(t-l)j tM- (3.11) 

The significance of this is that an identity of Sylvester (see |2S1) gives that if 

T3[0] = 1, (3.12) 
then the general solution of H3.10|l is given by 

f3H=det[<5^+'=f3[l]l . (3.13) 

Furthermore, restricting the parameter space so that 02 = (which corresponds 
to a chamber wall or reflection hyperplane in the afhne root system), it 
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has been shown by Okamoto |SZI that f3[l] is given in terms of a solution of 
the Gauss hypergeometric equation. Using integral solutions of the latter, the 
formula H3.13I) was taken as the starting point by Forrester and Witte in an 
extensive study of multidimensional integral forms of the r-function sequence 
fafn]. In particular, results relating to averages of the form (|2.26(l . equivalent to 
the following were established. 

Proposition 3. Define 

( TT(1 - ex[l0 ( 2 sin ^) ( —1^ " (1 - ue^'' f^^) (3.14) 
where zi = e*^' , < 0; < 2ti , and 

\0, 6iiJ. 

Let 

h=(^^/^{N+u-^i),uJ + ^/2{N+u+n),^/2{N-uo+^x),-^x-^/^{N+uo+^i)^, (3.15) 
and write 

Ci = e'2 [b] + tiN, C2 = 1/262 [b] + ^lN 

(recall the definition o/e2[b] and 62 [b] from Proposition 1). The Pyi system with 
parameters \3.15]) permits the r-function sequence 

T3[Ar] cxu^^/2AAr(u;c^,Ai;0 (3-16) 
where the proportionality factor is independent of u and furthermore 

Ciu ~ C2 + u[u - I) log Ai^ In; uj, ^j,; = hyUu; b) (3.17) 
au 

where /ivi(t;b) is an auxiliary Hamiltonian 1^,^.6]} for the Pyi system with pa- 
rameters f3.15\) . Consequently {3.17}) satisfies the Pyi equation in a-form 7| ) 
with parameters ^3.15}) . 



To relate to we note that 

)(^ 



(1 - 2x[l0 ^ (l-ue^'0L=e-^* = 2|sin. 



Consequently 

P%+ii^;0) = ^A^{e'"-/^;y2,y,;2) (3.18) 
where we have used the fact that p'^^^{x; 0) is even in x. 
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The choice of the parameters in H3.18|l corresponding to p'^^i imphes a special 
structure to the r-function sequence (pnH|l . First substituting (^1^ in l|XT^ 
shows we are considering the Pyi system with parameters 

b = ( VsTV, 1 + V27V, - 1 - i/^N) . (3.19) 

As noted above, tz[1\ satisfies the Gauss hypergeometric differential equation. 
The parameters in the latter are related to the parameters b by 

a = &i + 64, = 1 + 63 + 64, c = 1 + 62 + 64- 

Substituting the special values H3.19|l we see that in particular c = 1, which is 
the condition for the existence of a logarithmic solution at the origin [u — Q). 
For general TV, t^[N] then corresponds to a generalisation of this logarithmic 
solution of the Gauss hypergeometric equation. To illustrate this point, we note 
that with b given by H3.19|) . according to (|3.1t)|) and H3.18|l we have 

r3[7V](u)cx^/^/V^+i(x;0) 

Recalling H2.17|l and H2.18(l we see 

T3[l](u) cx {u + l)v + 2{u-\) 

T3[2](u) cx 4(1*2 ^^^^ ^^^2 ^ ^2(u - l)(u + 1)?; 

- u-^{u - 1)2(^2 _ I4u+ 1) 

where w = ttz — log u which exhibits the further structure of being a polynomial of 
degree N 'm v, and a Laurent polynomial in u of positive degree N and negative 
degree iV — 1. 

The Pyi system with parameters (I3.19f) also permits a r-function sequence 
which is strictly a polynomial. To anticipate this we relate (|3.14() to the free 
Fermi average lOnjl by noting 



_1 \ 1/2 



l-we'^OU^e— * =2sin- 



and so deducing 



{x; 0) = iA^(e2"-/^; 1/2, 1/2; 0). (3.20) 



Recalling (|3.16fl and H2.31|l we see that this corresponds to the r-function se- 
quence 

N-l 

T3[N]{u) CX ^ W^'. 

3=0 

This class of polynomial solutions is a special case of the generalised Jacobi 
polynomial solutions identified in |35| . As a final remark on the theme of special 
classes of solutions to the Pyi system, we note that the specification of the 
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parameters H3.19II is a particular example which permits elliptic solutions jSEHSl 
I17| . More generally the latter occur when 

ii = l + 63-64 = 2 + AreZ 
t2^bi + b2 = l + N eZ A 

<3 = &i - 62 = -1 e Z ^ 
^4 = 1 + 63 + 64 = e Z 



Substituting H3.18|l into H3.17|l of Proposition 3 and replacing iV by — 1 
throughout gives the sought evaluation of p%{x; 0) in terms of a solution of the 
Pvi equation in cr-form. 

Corollary 1. Define 

crjv(u) := u{u - 1)^ logp^(x; 0)\^2^^./L^^ (3.21) 

so that 

(rx/L J, \ 

The quantity (Jn{u) satisfies the particular Pyi cr-form differential equation 

+ [gn ~{u~ l)a'^ + l]{4f7^((7Ar - ua'^) - {N^ - l)[aN ~ {u - l)a^v]} = 

(3.23) 

subject to the boundary condition 

^n{u) ^(w-l)' + ^ 7^ i{u-lf + ... (3.24) 

The formula 1^3.2'J\) also holds for p^^^ [x; 0), except that a^f {u) is now subject 
to the boundary condition 

Proof - This is immediate from Proposition 3 and (|3.18() . H3.20|l . except for 
the boundary conditions. The latter in the free Fermi case follows by substi- 
tuting the exact evaluation (|2.31() in H3.21|l . Use is also made of the free Fermi 
density matrix exact evaluation H2.31|l to deduce the boundary condition in the 
impenetrable Bose gas case. Thus according to (|2.36|) - (I2.38() we have 



p^(x;0) = p^^^(a;;0)-2 / det 



dx-). 



(7V-i)(iV + l)/7r.T\2 {N -1)N{N +1) /TTx\3 



xZqP^V 6 \LJ ' 97r \L 



(3.25) 
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where the second hne foUows after substituting H2.31(l and expanding the first 
term to O(x^) (this term only contains even powers of x), and the second term 
to its leading order, 0{x^). Finally we substitute (|3.25|l in (|3.21|l to deduce the 
expansion (|3.24|l . □ 
One immediate consequence of Corollary 1 is that it allows the small x ex- 
pansion to easily be extended. Thus it follows that the corrections to H3.25|l at 
order and are 

(N - 1){N + 1)[3N^ ~ 7] /7rx\4 (iV - l)iV(7V + l)[lliV2 - 29] /7rx\5 

^ 360 vTJ 13507r vT/ 

(3.26) 

The results (|1.5|l . (|1.6() of Jimbo et al follow simply from our results 
(E!23)- Thus defining av{t) = HmAr^oo crAr(e2**/^) we obtain 1T3| from 
while substituting u ^ e^**/^ in replacing aNie'^'*^^) with cry(t) 

and equating the leading order terms in N (which are 0(1)) to zero gives (|1.6|) . 

The boundary condition l|1.7|l corresponds to the scaled limit of H3.24|l . 



3.2. Orthogonal polynomials on the unit circle. A feature of a number of re- 
cent studies |1UII11|[T1I5IB) relating Hankel and Toeplitz determinants to Painleve 
transcendents has been the characterisation of the former not only as the so- 
lution of nonlinear differential equations, but also as the solution of nonlinear 
difference equations. Here we will show a difference equation characterisation is 
also possible for p%{x; 0). 

For this purpose we adopt an orthogonal polynomial approach, similar to that 
used in |19j . The characterisation of the density matrix as a Toeplitz determinant 
with a non-negative and bounded symbol (|2.13I2.14|I immediately implies an 
underlying orthogonal polynomial system defined on the unit circle. The weight 
appearing in (|2.14(l is the special case a = b = 1/2 of the generalised Jacobi 
weight 

u;(z) = ^|l-|-^|2°|l^-uz|2^ a,beC,zeT (3.27) 



2tt 

where C is the normalisation 

C f dz 
2^ 



— |l + zp°|l + uz|2^ = l. (3.28) 

T 



Associated with (|3.27|l is a system of orthonormal polynomials {0n(2)}^^o,i,...' 



— W(z)0„(z)(/)„(z) = Sm.n- 

Y tz 



In obtaining a recurrence relation for 

dot 

and thus since 



—w{z)z''-^ 

T 



j,k=l,...,N 



P%+i{x; 0) = -l^DM-i\a=b=i/2 (3.29) 
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for the density matrix, one focuses attention on the leading two coefficients k„, 
In in (j)n{z), and the traihng coefficient (/'n(O), 

0„(z) = K„z" + Z„z"-i + . . . + 0„(O). (3.30) 

The relevance of k„, (/'„(0) are seen from the Szego relations 0J 



2 _ ^n-l 



fc— 



which show in particular that 

l-\rN\'^ ^''r^" , ,„:=nil-. (3.31) 



l2 _ Dn-2 Dn ^ _ 0n(O) 



'N~l 



We will see that for the weight H3.27|l , the Freud equations - which are recurrence 
relations among the successive coefficients k„,(/)„(0) ~ have a special structure 
which leads to a recurrence equation for r„, and thus according to H3.29(l and 
(lOTTl . for /9^(x;0). 

Proposition 4. Consider the special case a — b — 1/2 of {3.27)) , in which case 
according to \y.2y\) the relation 1^3.8 1}) reads 



1-|.^P = ^4±H^. (3.32) 



The ratios rn, and thus via 1^3. 32\) the successive density matrices, are determined 
by the third order difference equation with respect to N 

2 cos ^ + 2fN+ifN ^ ^7''^+' [{N + 3)rN+2 + {N + 1)Pn] 

^ ^^+1 , (3.33) 

l-fi 



rN 



^ [{N + 2)fN+i + NfN- 



1 



where rn := e^^'^^ ^/-^)"f„ G M. The initial members of this sequence of rn 
required to start the recurrence are 



ro = 1 (3.34) 
TT — 2ttx/L + sin(27ra::/L) 
^(tt — 2ttx/L) cob{ttx/ L) + sin(7rx/-L) 



(substituting these values in yj. 33\) with N ~ allows r2 to be computed). Also, 
the initial members of the sequence are specified by \2.1b]) and \2.n\j (with 
these values given pg is computed from 1^3. with N ~ 1). 

Proof - Magnus has found a recurrence relation ^2 ^oi the ratios r„ applicable 
to the generalised Jacobi weight (|3.27() 

/ , -, , , ,x , / -, , , ,x- uln/Kn + ln/Kn-n{u+l) 

{n + l + a + 6)r„+i + (n-l + a + 6)ur„_i = j — (3.36) 

1 rn. 
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However this involves both r„ and Z„ and we require a further relation to deter- 
mine In- This relation is 



In _ L (a - b)n 

u— = — — — — (w- 1), 

Kn Kn n + a + 



and follows from telescoping the identity 



(n + a + b+l) 



In+l 



n+1 



which in turn is derived by evaluating 
dz 



-{n + a + b) 



(3.37) 



= (a-5)(S-l), 



(1 + z)(l + uz)w'(z)0„(z)0„(z) 

T «2 

in two different ways. Then H3.33|l follows after setting a — b — y2 and extracting 

a phase factor of e'^^Ci-^/- ^)" 

The initial conditions can be determined by a Gram-Schmidt 

type construction of the orthonormal polynomials {(t>n{z)}n=a.i,...- First, due to 
the normalisation of the weight H3.28|l we have (/'o(z) = 1 and thus H3.34(l follows. 
With 

f dz 

if, 9) := / —w{z) f{z)g{z), 

Jf IZ a=b=l/2 

the orthogonality (z), (/)o(2;)) = 0, explicit value 4iq{z) — 1 and H3.30() give 
Ki(z, l)-|-(/)i(0) — 0. The value of (z, 1) can be read off from (|2.15|l . thus implying 

(ESS)- □ 

We note that the special cases corresponding to x at either the end points 
[x — and thus it = 1) or the midpoint {x = L/2, and thus u = —1) allow 
simple explicit formulas for the r„. Thus we have jl9| 



(-1) 



N 



1 



iV+ 1 



-1, rN=2p = 



1 



(3.38) 



N + 1' 



rN=2p+i = 



which clearly satisfy (|3.33l) . The density matrix at these points also has a closed 
form evaluation, 



p^(0;0) = 



N 

T 

i f^Y'\(.-i)(2r,-i) G{p + 2)G^iP+l)Gip) 
L\tt) G^{2p+l) 

i (iY'\pi2p-i)GHp±})GHl±i)_ 
L\tt) G2(2p + 2) 



where G{x), the Barnes G-function, has the explicit form G{x) — (x — 2)!(2; — 
3)! ... 1! for X G Z>2. Here the former evaluation follows from the general fact 
that at coincident points the density matrix is equal to the particle density, while 
the latter makes use of results from IW\ on the explicit form of the k„'s for the 
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weight (jS2ZI), = K2n+i = (2n + l)!/2^ "(n!)^ \/n + 1, in the case u = -1, 
a = b. The small x expansion of p'^{x;0), H3.25|l and H3.26|l . substituted into 
(I3.32|) allow the corresponding small x expansion of f„ to be computed up to a 
sign, which in turn can be determined using (|3.38|l . This shows 

A feature of the explicit forms (|3.38|) is that |rjv| ^ as iV — > oo. According 
to (|3.31|l this is a necessary condition for the convergence of p%{x; 0) as iV ^ oo. 
Here we note that with |rjv| small the difference equation (|3.33|) simplifies to read 



2cos^ = S^+l-i3^., j,^.^^iN+^}rN^l±Nr^^ 
L rN 
It follows from H3.39|l that Bm ~ 2N costtx/L + C, where C is independent of 
N. Noting that this implies ^ 2{N + 1) cosnx/L for N large, we thus obtain 
the recurrence 

{N + 2)Pn+i +NrN-i = {N + l)fN2cosTTx/L, (3.40) 

which with An '■= (N + l)rN, cosirx/L := t reads 

AN+i+AN-i=2tAN. (3.41) 

This is precisely the three term recurrence satisfied by the Chebyshev polyno- 
mials HJI). Thus H3.33|l can be regarded as a non- linear generalisation of the 
Chebyshev recurrence (|3.41() . 

Although it is not obvious from the derivation, the equations ()3.32(l and H3.33|l 
remain valid in the free Fermi case. This can be seen by substituting the exact 
evaluation l|2.31|l into H3.32|l to deduce 

sinTTx/L 



sin7r(ri + l)x/L 

and then verifying that this is an exact solution of H3.33|l . Unlike f„ in the Bose 
case, (|3.2ll does not obey the inequality |f„| < 1 for all x. 

As our final point of the difference equation, we remark that recently Adler 
and van Moerbeke 1 have constructed essentially the same pair of coupled re- 
currences ()3.36(l . H3.37|l from their theory of the Toeplitz lattice and its Virasoro 
algebra. In the particular case at hand their weight is specialised to a = /? = 1, 
and ^"■^ = u = e^^^'-^/L, Their variables are related to ours by x„ = f„ and 
through the use of their relation (0.0.14), which is the analogue of (|3.37() . then 
one can show ?/„ = x„. The other recurrence in their work, (0.0.15) is the ana- 
logue of l|3.36|l and can be shown to lead to 

(1 - x^) [{n + 2)a;„+ia;„_i + n + 1] - (1 - x^_i) [{n - l)xnXn-2 + n] 

= 1 + 2Xn~lXn C0S{-KX/L) + x'^_ix'^ 

- 1 + (1 - xl){3x2 + 2) - xi{xi +2cos{ttx/L)). 

Now by clearing denominators and rearranging (|3.33|) one can recover the first 
five terms of the above relation. Furthermore by using the initial conditions of 
the recurrence (|3.34|l . H3.35|l one can show that the sum of the last three terms 
is identically zero and thus the two forms are the same. 
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4. Painleve-type evaluations of pn{l{x);x) 

4-1. Fredholm formulation. While the unitary average (|2.26|l defining p^j^-^[x\ 0) 
is a known r-function in the Painleve theory, the same is not true of the averages 
(12.27(1 - (|2.29|l . Indeed the density matrices in these cases are genuinely functions 
of both X and y. These variables play the role of time in the Hamiltonian formu- 
lations of the Painleve equations, so there being more than one time variable, 
we are taken outside this class. However, with y — l{x), where l{x) denotes the 
reflection of x about the centre of the system (thus i{x) = —x for the harmonic 
well case, and l{x) — L — x for the case of the Dirichlet and Neumann boundary 
conditions) we again have a function of one variable. Although this cannot be 
recognised as a single r-function, it turns out that we can formulate the cal- 
culation of pN{i{x);x) so that it is expressed in terms of quantities known in 
terms of Painleve transcendents from random matrix theory. For this one makes 
use of a classical operator theoretic interpretation of H2.l-i8|l relating to Fredholm 
integral equations ,30^ . 

It is the latter formulation which has been used in the pioneering work of 
Jimbo et al. |22| on the evaluation of the bulk density matrix in terms of a 
Painleve V transcendent, and the generalisation of this result by Its, Korepin 
and coworkers |20II28| to the temperature dependent bulk density matrix. The 
key point is that with Kj denoting the integral operator on J = [x,y\ with 
kernel (|2.37|) . and R{a,b;^) denoting the kernel of the resolvent operator R := 
^Kj{l - ^Kjy^, it is true in general that (see e.g. pUll^ ) 



^[x.,] ( 5 ; = det(l - ^Kj)R{a, 6; (4.1) 

(the quantity '^ix,y] is called the first Fredholm minor). Now in the harmonic well 
case and the cases of Dirichlet and Neumann boundary conditions (the latter 
two after the change of variables cosirXj/L k-* xj) the wave function is of the 
form with 

2f \_ J e^^ , harmonic well (ao\ 

^ 1(1- a;2)±i/2^ Dirichlet and Neumann ^ ' ^ 

These weights have the property of being even is x. This implies a special struc- 
ture to 1(4.111 if J is also chosen to be symmetrical about the origin, J = [—x, x] 
say. Thus a consequence of g^{x) being even is that the orthogonal polynomi- 
als Pj{x) are even for j even and odd for j odd, and this from (|2.37() implies 
K(a, b) = K{—a, —b). Using this latter property, and with J — [^x, x], it is true 
in general that (see e.g. |32]) 

logdet(l - Ci^j) - -2i?(x, x; e). 

ax 

Antidifferentiating and substituting in 1(4. l|l with [x,y] [— Xjx], then substi- 
tuting the result in (|2.38|) shows 



Pn{-x;x) = R{-x,x-^)exp(^-2 J R{t,t;^)dt^ ^ ^. (4.3) 
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The crucial point of the formula 1)4.3(1 is that the quantities R(—x, x) and 
R{t,t), for Christoffel-Darboux kernels corresponding to the weights 1(4.2(1 have 
previously been calculated in terms of Painleve transcendents as part of studies 
into gap probabilities (interval J free of eigenvalues) for random matrix ensem- 
bles, the GUE in the case of the harmonic well, and the JUE with a = 6 = ±l/2 
in the case of Dirichlet and Neumann boundary conditions. Although 14.2(1 and 
((4.3(1 have general validity, the specific integrable nature of the kernel 1(2.37(1 20 
is essential for this characterisation. In the latter case the quantities in ((4.3(1 
were studied in j49[ . but with J — (—1, —a;] U [x, 1) rather than J ~ [—a;, x]. To 
overcome this difference in detail, we note we can rewrite ((2.33(1 to read 



c 



dxN 9^{xn) 



N 

X JJ(a; - x/)(j/ - a;;) {xk - xj)'^ 

1=2 2<j<k<N 



Repeating the working which led to 1(4.3(1 then shows 

Pn{-x; x) = i-lf-'Ri^x, X- C) exp ( - 2 / R{t, t; dt) | (4.4) 



where R now denotes the kernel of the resolvent operator R — j{l ~ £,Kj) ^, 
J :— (— oo, — x] U [x, oo). 



4.2. Evaluation of p^{~x;x) and {L — x; x) . Let us begin by specifying the 
quantities in 1(4.3(1 in the harmonic well case. From the above discussion, this 
corresponds to the interval J being eigenvalue free in the GUE. For this matrix 
ensemble, the gap probability for both the interval J ~ [—x, x], and the interval 
J = (—00, —x] U [x, 00) have been studied [43II49[ . allowing us to use either ((4.3(1 
or ((4.4() to deduce an exact expression for p^{—x; x). For purposes of specifying 
the boundary condition in the corresponding differential equation, it is most 
convenient to use the latter. 

Proposition 5. For the impenetrable Bose gas on a line, confined by a harmonic 
well, we have 



(-x; x) = E(a;; exp ( - 2 / R{t; ^) dt 



5=2 



where with h a/s^ — 2R' , R satisfies the equatic 



sR" + 2R' ^ 2s{s -h)- 2h^{R + sR'f - As'^{s - h)R - 2Ns^{s - hf 
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while R satisfies the equation 



sR" + 2R' + 8NsR - 24s^R^ 



4(s - 2Rf(^{R + sR'f + ms'^R? - 165^7?^^ . 



The corresponding boundary conditions are 

2N-1 



N-1 



7ri/2(iV_ 1)1 



□ 



The two resolvent kernels are not independent being related by 
d 



ds 



i?(s;C) = -2si?(s;0-2i?^(s;C) 



so that h = s + 2R{s; ^). Both resolvent kernels have been reduced to a particular 
Pv transcendent w{x) with parameters 



V2, 



where e = ±1. The reductions were found to be 



R = 



1 



%^W{W - 1)2 



2x—w + N{w - if + {2x - l)w + 1 
dx 



X < 2a;— —ui — Nlw — if - (2a; + l)w + l 
dx 



R = - 



2w{w - 1) 



t—w — w 
dx 



N{w-l)+e 



where s"^ — x. 

In the case of Dirichlet and Neumann boundary conditions we require R{—t, t) 
and R{t,t) for the symmetric JUE with {—1,—t] U [t,l) eigenvalue free. The 
differential equation satisfied by R{t,t) is known from 35!, but the equation for 
R{—t,t) was not made explicit in that work. We therefore give some details of 
the required calculation below. 

Proposition 6. The impenetrable Bose gas on the finite interval [0, L] subject 
to Dirichlet or Neumann boundary conditions at the ends has the density matrix 



Pn [L-x;x} 



\ sin ^-^R^,^"" is; C) exp (^-2 j' dtR^^'^it; 
where a{s) :— (1 — s'^)R^'^ (s; £,) satisfies 



F 



[{N + afs + 2s(t' - (1 - s^)a"] + (1 + s^)F + 2{N + a)s 



- {2(1 - s^)a - 2{N + a)s^[{N + a)s + F]- s[{N + a)s + Ff) 

-As'^[{N + a)s + Ff {N{N + 2a) + 2s<t + 2{N + a)s[{N + a)s + F]} 

(4.5) 
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where a — ±1/2 and F :— ^ (N + (xy-s^ — 2(1 — s^)cr' and with the boundary 
condition 



-M-ir-'pTi^;^) (4.6) 



^ (£ - 1)^-^ iN + a)r{N + 2a + l) _ 
s-i^^ ^ 2"+i(iV- l)!r(a + l)r(a + 2)^ ^ 

and i?o ■— R^'^{s',0 satisfies 

- s^)R'^ + 2(1 - s2)(l - 2s^)R'o 

+ 8sRo{2sRo - {N + a)/2)(2si?o - iV - a) - 2(1 - + 2a2)si?o| 
= A{-{N + a)s + 2{l + s^)Ro} 
X {(1 - s2)2(i^„ + si?^)2 + 4(si?o)2[(2si?o - iV - a)2 - a^lj (4.8) 

with the boundary condition 

Rn^-.o ^(e-i)^-^i^r (4.9) 

1 r(N + 2q: + 1) , , , 

Proof - Both cases come under the symmetric Jacobi weight discussed in 
H^j however this study doesn't contain sufficient details for our purposes. From 
Proposition 8 of the above we can derive an integral of the system of differential 
equations. Combining (4.26) and (4.27) we have 

(1 - s^)iqpy - [f3o + u(2ai-l)]/ + [70 - w{2ai + l)]q^ = (4.11) 

and we can rewrite (4.23) as 

a~[(3o+u(2ai-l)]p^~[-fo~w{2ai+l)]q^-2aisqp-2s-\l~s'^)q^p^ = 0. (4.12) 

Adding to 2ai times (|TTT|l and employing (4.24,25,28) we find the sum 

to be a perfect derivative, and given that Ro,a,u,w all vanish as s 1 the 
integral is 

[(3o + u(2ai-l)][7o - w(2ai + l)] - 2sa - 4ai(l - s^)qp = /3o7o ^ N{N + 2a). 

Given this integral we can follow the procedure in the proof of Proposition 1 1 
in 48 to derive second-order differential equations for a and Rq, and equations 
(14.51) . 1)4. 8|l are the results. The boundary conditions arc found from the kernels 

NiriN + 2a + l) 

KN{t,t)^ 22"+ir2(7v + a) ^^-'^ 



X 



K^ft t) ( )^-' N\riN + 2a+l) (l-t^)° Ko) 

) 2^-+^r{N + a)r{N + a + l) t ^^-^ ^*>^^ ^ 
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where the P^'^\t) is the standard Jacobi polynomiaL 
Again the two resolvent kernels are connected by 



□ 



—a = -2{N + a)sRQ - 2(1 - s^)RI 
This system was solved in terms of Pyi transcendent w{x) with parameters 



V8, /3 = -V2(iV + a-e/2)2, 7 = 0, 5 = 1/2(1 -a^) 



according to 
i?o(s) 



-w 



dx 4y^(x — 1) 



2{x^l)-^w-{w-l){w + x) 



[e{w + x)~ 2{N + a)x] 
2 



8y/xw(w — l)(u' — x) 
^/x{w - 1) 



2'w{w — x) 



[N{N + 2a)w — (N + a)'^x] 



where s = ^Jx and e = ±1. 

We remark that the thermodynamic form of (x\ y) for x and y fixed but 
general has been studied in the spirit of the paper of Jimbo et al. by Kojima 
\27\ . although the final characterisations obtained (in terms of an integrable 
differential system) does not appear to be amenable to computation. 



4-3. Thermodynamic limit. In the thermodynamic limit, Jimbo et al |22j were 
able to identify the first Fredholm minor directly as a r-function, and so had 
no need for the formula (|4.3|) . Nonetheless it can be used to evaluate Poo{x; 0) 
in terms of a solution of (|1.8|) . as we will now demonstrate. First, repeating 
the working which lead to (|4.3|l shows that with Kj the integral operator on 
J — {—x,x) with kernel K{x,y) = sm{x — y)/Tr{x — y) and R — ^Kj{l — ^Kj)^^, 



Pooi-x;x) = TrpQRooi-TTpQX,Trpox;(,)exp 



Roo{t,t;^)dt 



«=2 



Furthermore, we know from 22. that 
t 



:i?oo(V4 t, 1/4 t; = hy{-it- (0, 0, 0, 0)), 



(4.13) 



where /iv(i;v) satisfies (|1.8|l . subject to the boundary condition 

hy{-zt;{0,0,0,0)) -^^ - 

t^o 2tt 47r^ 



and we know too that 
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(see e.g. Consequently Poc{x; 0) can be expressed in terms of the transcen- 

dent l|4.1l-{(l according to 



Poo (a;; 0) = vrpo 



d hv{-it; (0,0,0,0))V/^ 



dt 



X exp 



/iv(-m; (0,0,0,0)) 



du 



t = 27rpQ X 

e=2 



(4.14) 



This is to be compared against the evaluation H1.5|) due to Jimbo et al |22| . 
with the substitution (|1.9() and the boundary condition l|1.7|l generalised to be 
consistent with H2.38|l . 



Poo(a;;0) = poexp 



* -t 



du 



5 = 2 



(4.15) 



/iv(-*i;(V2,-V2,V2,-V2)) 



1 



t^o 2 12 487r 

Equating the logarithmic derivatives of (|4.14|) and H4.16|l gives the identity 



(4.16) 



hvis; (V2,-V2,V2,"V2)) 



s d 



'h + hvis; (0,0,0,0)) + --log 
2 ds 



d /iv(s;(0,0,0,0)) 



It remains to derive this directly from the Painleve theory. In this regard we note 
that an identity with similar characteristics for the Pn system can be deduced 
from a classical result of Gambier, while the analogous result for the Pm system 
has only recently been found |47|. 



5. Discussion 

We will conclude with a discussion of our results. In relation to p^(a:;0) the 
recurrence relation H3.33|l allows rapid and stable tabulation for very large values 
of N for all x E [0,L). Hence numerical evaluations of the Xk in 1)1. 3|l can be 
carried out. Although for fixed k the leading behaviour in A'^ of are known 
from 1)1.4(1 . of interest is the convergence of the Xk (appropriately scaled) to their 
thermodynamic value. Furthermore, the differential equation characterisation 
of p%{x;0) given in Corollary 1 is well suited to generating the power series 
expansion about a; = 0. This relates to the behaviour of Xk for k large. A 
comprehensive study of such issues will be discussed in a forthcoming publication 

Our results show a remarkable Fermi-Bose correspondence at a mathematical 
level of characterising the density. It goes back to Girardeau that up to the sign 
under permutation of the coordinates, the ground state wave function of the 
Id free Fermi and impenetrable Bose systems are identical. This means that 
quantities depending only on the absolute value squared of tpo are the same 
for both systems. The density matrix is not of this type, and so distinguishes 
the two systems. Nonetheless, we find that the same differential and difference 
equations characterise p%{x; 0) for both the Bose and Fermi systems - they are 
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only distinguished at this mathematical level by the boundary conditions. For 
Pqo{x; 0) this same property was already observed by Jimbo et al 1221 . One would 
like to be able to use the differential equation to see how the different prescribed 
small X behaviours imply different behaviours at large x. In the case of Poo{x] 0), 
Jimbo et al |22 were able to do this and so obtain the expansion 



which was first derived by Vaidya and Tracy |44II45| . A challenging problem is to 
use Corollary 1 to deduce higher order terms in the expansion H1.4|l of p%{x; 0) 
for large N, with x/N fixed. Note that in the Fermi case this expansion can be 
read off from H2.31I) . 

For the density matrix p^{x;y) (and similarly p^'^{x;y)) we have the pre- 
sented determinant formulation specified by (|2.19|) . H2.22|l . as well as the formu- 
lation H2.27|l as an average over the eigenvalues of the GUE. Both these forms 
are suitable for the numerical computation of p^{x;y) for general x and y. In 
the special case that (x; y) i— *■ (— x; x) we have given an explicit functional form 
in terms of the transcendents related to Py • This provides a much more efhcient 
numerical scheme for the computation of p^{—x; x), and will provide a valuable 
test for the accuracy of Monte Carlo evaluation via H2.27(l . 

With this achieved, the next step in the determination of the occupation 
numbers is the numerical solution of the integral equation (|1.2|) . In addition to 
the numerical evaluation of the occupations, one would like to determine the 
exact leading asymptotic behaviour of Ai, i fixed, from an asymptotic expression 
for p^{x;y) analogous to l|1.4(l . Thus we seek the asymptotic expansion of the 
determinant specified in (|2.19l) . (|2.22|) . which can be considered as a Hankel 
generalisation of a Fisher-Hartwig type Toeplitz determinant, or equivalently 
the asymptotic expansion of the random matrix average l(T?7j) . On this latter 
viewpoint of the problem, we recall that the original Szego theorem for the 
asymptotic form of the Toeplitz determinants with smooth, positive symbols 
has been proven by Johansson 2$ starting from the analogue of the random 
matrix average (|2.26|l and then generalised to an analogous theorem for Hankel 
determinants relating to Jacobi averages (24j . We note too that some mappings of 
Fisher-Hartwig symbols in the Toeplitz case to analogous symbols in the Hankel 
case are known 

Another aspect of the present work which provides the beginning for future 
studies is our derivation of the recurrence relation H3.33|l using orthogonal poly- 
nomial theory. As we have commented, the recurrence obtained via this method 
coincides with recurrence obtained by Adler and van Moerbeke using methods 
from soliton theory. Of course one would like to understand the underlying rea- 
son for this coincidence. Also, recent works of Borodin and Boyarchenko 0,0], 
starting from a formulation in terms of the discrete Riemann-Hilbert problem, 
provides alternative recurrences involving discrete Painleve equations for closely 
related Toeplitz determinants. It remains to understand the relationship between 
the different approaches. 

Note added: Subsequent to the completion of this work an asymptotic form 
analogous to 1)1.4(1 has been derived in [2j for the harmonic well case, and from 
this it is deduced that, as with periodic boundary conditions, the for fixed k 
are proportional to ViV. 




(5.1) 
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